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Abstract 

The interface tension between Z(N) vacua in a hot SU(N) gauge 
theory (without dynamical fermions) is computed at next to leading order 
in weak coupling. The Z(N) interface tension is related to the instanton 
of an effective action, which includes both classical and quantum terms; 
a general technique for treating consistently the saddle points of such 
effective actions is developed. Loop integrals which arise in the calculation 
are evaluated by means of zeta function techniques. As a byproduct, up 
to two loop order we find that the stable vacuum is always equivalent to 
the trivial one, and so respects charge conjugation symmetry. 



I. Introduction 

In the absence of dynamical fermions, SU(N) gauge theories possess a global 
Z(N) symmetry associated with the center of the gauge group [Q. Confinement 
implies that at zero temperature the vacuum is Z(N) symmetric, but at nonzero 
temperature there is a phase transition to a deconfining phase, where the Z(N) 
symmetry is spontaneously broken 0. In the deconfined phase, a system of infinite 
extent falls into one of the N degenerate vacua, but in a finite volume bubbles of 
different vacua form || , separated by domain walls. The action of these domain walls 
is proportional to the interface tension, and so controls the dynamics of large Z(N) 
bubbles. 

We previously discussed how to compute the interface tension a in terms of the 
temperature T and the coupling constant g Q. In weak coupling, 

a ^-jy Zfl- (15.27853...)f^V (1.1) 
3V3N g V 167r J 

The term at leading order, T 3 /g, is the result of Q; the principal result of this paper 
is the computation of a at next to leading order, gT 3 . In this expression the coupling 
g represents a running coupling constant at a temperature T; the exact relationship 
of this g to the bare coupling constant is given in ( |4.11| ). 

The physical applications of our results is at best indirect, for it only applies to 
a world devoid of dynamical quarks. It can be compared to numerical simulations of 
lattice gauge theories ||; indeed, it was these measurements of a which led us to ask 
if it is computable in weak coupling. 

Nevertheless, the methods and techniques which we have developed to compute the 
interface tension are, we believe, of general interest. The calculation of the interface 
tension reduces to an instanton problem in an effective theory. This effective theory 
is one dimensional P, 0, as it describes the profile of the interface in the direction 
perpendicular to the domain wall. Remember how a standard instanton calculation 
proceeds ||: the instanton is the solution of the classical equations of motion, with 
the action for the domain wall proportional to 1/g 2 . The corrections to the classical 
action are given by expanding the classical action in a background instanton field. 
Corrections at next to leading order are given by integrating over terms of quadratic 
order, which gives corrections to the action of order one. 

Contrast this with the interface tension in ( |1 . 1|) . The action of the domain wall is 
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a times the transverse volume. By its mass dimensionality, a is proportional to T 3 , 
while the constants, and the factors of N, result from the detailed path the Z(N) 
interface takes in the space of SU(N) gauge fields. What is surprising is that a 
starts out as 1/g in weak coupling, for if the Z(N) interface were a solution to the 
classical equations of motion, it would start instead as 1/g 2 . This is because the 
effective action (in one dimension) which controls the Z(N) interface is the sum of 
the classical action plus a quantum term, obtained by by integrating out fluctuations 
at one loop order. In the effective action the classical piece acts like a kinetic term, 
and the quantum piece like a potential. The Z(N) instanton is a stationary point 
only of the full effective action, and this balance between classical and quantum terms 
transforms the usual factor of 1/g 2 into a 1/g. 

It is then unclear how to compute corrections to a beyond leading order. Surely 
one cannot compute blindly: if the effective action is expanded in fluctuations about 
a background instanton field, and these fluctuations integrated out, how is double 
counting avoided? That is, how to differentiate between the quantum fluctuations 
which generate the effective action in the first place, from the quantum fluctuations 
which are properly included in the expansion about the instanton? 

In this paper we solve this problem for the Z(N) interface in a manner applicable 
to arbitrary effective actions. The method generalizes what is known as the "con- 
strained" effective potential |7j]. Here we use it to reduce the four dimensional gauge 
theory to an effective scalar theory in one dimension. The method is trivial in design: 
a delta function constraint for new degrees of freedom is inserted into the functional 
integral. The original degrees of freedom are then integrated out, producing an effec- 
tive theory for the new field. There is obviously no problem with double counting, 
since extra degrees of freedom are introduced in the first place. 

For the interface tension, this method shows that the leading corrections are of 
order g 2 times that at leading order. These effects are due entirely to corrections in 
going from the four dimensional theory to the effective, one dimensional theory. They 
enter both for the kinetic and potential terms in the effective action. For example, 
the corrections to the kinetic term transform the bare into a renormalized coupling 
constant. 

The order parameter which distinguishes different Z(N) vacua is the Wilson line 
at nonzero temperature @. It is convienient to parametrize a vacuum expectation 
value for the Wilson line by giving the gauge field Aq a nonzero value; the effective 
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action for a Z(N) instanton is then related to the free energy in such a background 
field. The Z(N) instanton is slowly varying in space, so that previous results in a 
constant A field [§ - [14| , especially those by Belyaev and Eletsky [[□]] , Enqvist and 
Kajantie |I2 |, and Belyaev [JH| , can be used. The technical problem of computing 
the free energy in a background Aq field is done most easily by using zeta function 
techniques |T5| ]. 

Much of the interest in considering a hot theory with A ^ concerns the pos- 
sibility of the vacuum spontaneously generating an expectation value for A , and so 
for the Wilson line [PH^]]- If this happens for N > 3, it implies that the vacuum at 
nonzero temperature spontaneously breaks charge conjugation symmetry. For theo- 
ries without dynamical fermions, following Belyaev [|14[] we do not find evidence for 
the spontaneous breaking of charge conjugation symmetry at two loop order: up to 
the usual Z(N) rotations, the stable vacuum is the trivial state, with Aq = and the 
Wilson line equal to one. 

The outline of the paper is as follows. In sec. II we review the calculations of 
0]. Sec. Ill outlines the general calculation of the interface tension. The calculation 
of the interface tension at next to leading order in Feynman gauge is given in sec. 
IV. Sec. V considers the computation in arbitrary covariant gauges. There are two 
appendices. In appendix A we prove that the path chosen for the Z(N) instanton has 
minimal action for N = 3 and oo. Appendix B summarizes various integrals needed 
in a constant background A field. 



II. Z(N) interface at leading order 

We begin by rederiving the results for the Z(N) interface tension at leading order 
0| . This is not mere repetition, for here we use a more natural basis for the generators 
of SU(N) matrices than before. The subtleties of how to derive the effective action 
are deferred until sec. III. 

We work in euclidean spacetime at a temperature T, so the euclidean time r 
varies from to /3 = 1/T. In the spatial directions the system is a long tube, of 
length L in the z direction, of length L t in the two remaining spatial directions, Xt, 
with L ^> L t ^> (3 . The volume in the directions transverse to z is V tr = f3L^. 

A Z(N) interface is constructed by assuming that the system is in one Z(N) phase 
at one end of the tube, z — 0, and in another Z(N) phase at the other end, for z = L. 
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This forces a Z(N) interface along the z direction, with the action of the interface 
equal to the interface tension, a, times the transverse volume, V tr . As a practical 
matter, while L ^> L t , in the end both are taken to infinity. 

The Z(N) symmetry is determined by the the trace of the Wilson line in the 
fundamental representation, 



tril(A) = — tr \ V exp (ig [ Ao(x)dr 
N \ \ Jo j 



(2.1) 



V refers to path ordering. Let A Q have the value, 

2%T 



A%(x) 



gcN 



qt N , 



(2.2) 



where tw is the diagonal matrix 



t 



N 



1jv-i 








1 



-(N-l) J ^2N(N-1) 
In the presence of this A Q field the Wilson line equals 



(2.3) 



tr Q{A C 



^2wiq/N 



-2wiq 



(2.4) 



The trivial vacuum is A = q = 0, with trQ = 1; Z(N) transforms of the trivial 



vacuum occur for q = j, with trQ = e 



2mj/N. 



j is an integer between 1 and (N — 1). 



Thus the simplest Z(N) interface is constructed by promoting the parameter q in 
( |2.2| ) to a function of z, satisfying q(0) = and q(L) = 1, so that trQ = 1 at z = 0, 



and tr Q 



2ni/N 



at z 



In order to proceed we need a useful parametrization for the remaining generators 



of SU(N) . The diagonal generators are chosen similarly to ( |2.3| ). For example, 



HAT-l) 



y/2(N - 1)(N - 2) 



( 



\ 



l N-2 



: \ 



-(N-2) 




(2.5) 



Altogether there are — 1 diagonal generators, with tj from i = 2 to N. In SU (2) 
t2 is proportional to the Pauli matrix 03, while in SU(3) ts is proportional to the 
Gell-Mann matrix As- 
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For the off-diagonal generators we follow Belyaev and Eletsky ]TTj] and use a ladder 
basis. For example, 

l\ 



b N,l 



NA 



1 

71 



i 

72 







: / 





1 



(2.6) 



(2.7) 
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All elements not indicated vanish. Each diagonal generator ^ has 2(i — 1) ladder 
generators associated with it: tfj, with j running from j = 1 to i — 1. For a ladder 
generator only one element is nonzero, 

1 



1 



(2.f 



^2 ' v ^ ymn ^ 

We are working in the fundamental representation, so the matrix indices m and n 
vary from 1 to N. 

These generators form an orthogonal set, with the diagonal generators normalized 

as 

tr{Ut v ) = ~<T\ (2.9) 
and the off diagonal generators as 



. 



(2.10) 



Observe that the metric in the ladder basis is diagonal in the i and j indices, but 
off-diagonal in the ± indices. 

The advantage of the ladder basis is the simplicity of the commutation relations. 
For the generator tpj, since from ( |2.3| ) it involves the unit matrix in the first N — 1 
components, 1jv-i, the only nontrivial commutator of is with its associated ladder 
generators, ■. This commutator is just a constant times the same ladder generator: 



[t;v, t* 1 



N,j\ 



±cNt%j . 



(2.11) 



This relation is familiar from SU(2), from where up to overall constants, t<i is a 3 , and 
tf 1 are the matrices a ± . 
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The Wilson line in the adjoint representation can be computed using ( |2 . 1 1| ) . This 
is defined as 



ryab 
1 L adj 



^ = iv^TI (t a rexp(ig J\a (x),]\ A , (2.12) 



where "a" and "6" refer to the (N 2 — 1) adjoint indices, and [A , ] denotes the adjoint 
operator, [A ,]X = [v4 ,X]. For the constant A field of ( |2.2| ), the only nontrivial 
elements of the adjoint Wilson line are those involving the ladder operators t%j, and 
so its trace is 

tr£l ad3 (A d ) = 1 - j^-j (1 - cos(2Ttq)) . (2.13) 

The adjoint Wilson line is unaffected by the Z(N) symmetry: for the Z(N) degenerate 
vacua, where q is an integer, trfl a dj = 1. 

The classical action is 

S d (A) = J^dr Jd 3 x 1 -tr(Gl) , (2.14) 

where G^ v = <9 M A, — dyA^ — iglA^, A u ] is the field strength tensor. For the interface 
problem the field q of ( |2.2| ) is assumed to be a function only of z, the length along 
the tube. For this field the classical action reduces to 



S'(^)=V t , lw {N-l) Jd z ^ . (2.15) 

Using the classical action, a solution to the equations of motion is q(z) = z/L. There is 
no true interface, since the action vanishes like 1/L as L — > oo. But this is misleading, 
for classically there is no sign of the Z(N) symmetry, as all values of q degenerate. 

This classical degeneracy is lifted by quantum effects ||. This is shown by calcu- 
lating the action in the presence of the background field in ( j2.2| ). For the time being 
we assume that q is independent of z, and concentrate on the g-dependent terms 
which lift the degeneracy in q. With A^ = A d + A q ™, in background field gauge ]16 
the gauge fixing and ghost terms are 

S^(A, V ) = J^dr J d 3 x (^tr(D d A^) 2 + fj^-D^D^nj , (2.16) 

where _D M = <9 M — ig\A^ ] is the covariant derivative in the adjoint representation, 
D d = <9 M — ig[A d , ], and r\ is the ghost field. 
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For a constant field q it is especially easy to expand the full action, S d + S 9 ? , in 
quadratic order in the fluctuations A qu , and then integrate them out: 

Sr(A d ) = \trln (-Dl5» v + (l - ~) D d D^j - tr In (-D d ) . (2.17) 

The subscript on Sf u indicates the quantum action at one loop order. The first 
term on the right hand side is from the integration over the gauge fields, the second 
from that over the ghosts. Note that because A cl is assumed to be independent of 
spacetime, there are no terms in the inverse gauge propagator proportional to G d . 

This quantum action is independent of the gauge fixing parameter, £. To see this, 
note that the derivative of Sf" 1 with respect of £ _1 is 

9ST(A") KJ +( l-Q D ;'T:\). (2.18) 



ar* - 2 v " " \- D i v ' (--d?,) 2 

That is, the derivative is —D d D d times the gauge propagator in the background field. 
Normally this propagator is difficult to compute because the covariant derivative 
doesn't commute with itself. For a constant background field, however, it does, and 
so the ordering of the D d, s is inconsequential. Then (|2.18|) reduces to 

*M = §MD. (2.19) 

Thus the variation of the quantum action with respect to £ is a constant independent 
of the background field, which can be dropped. 

Adopting Feynman gauge, £ = 1, the commutation relations of the ladder basis 
reduce the color trace in Sf u to an abelian problem. As D d is the adjoint covariant 
derivative, it is independent of the background field unless it acts upon the ladder 
matrices ^ ■: 

D d t%j = (d T 2-KTqi) t% 3 = Dl t% j . (2.20) 

With the euclidean four momentum equal to (k°,k), at nonzero temperature k° = 
27mT for integral n, and the covariant derivative becomes 

iD° ± — ► k± ee 2nT(n ± q) . (2.21) 

The sum over n implicit in the trace includes both positive and negative values, with 
the sum over k°_ equal to that for A/j_. Hence the quantum action reduces to 

ST(A d ) = 2(N - 1) tr In ((fc° ) 2 + k 2 ) , (2.22) 
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k 2 = k 2 . 

This result is typical of loop effects in a constant background field. For the degrees 
of freedom along the ladder operators t^j , the propagators are as in zero background 
field, except that k° is shifted by a constant amount, to k±. The propagators for the 
remaining degrees of freedom are unaffected by the background field. 

From fl2~22l) , at one loop order the g-dependence of the free energy reduces to 
(N — 1) copies of that for SU(2). To isolate the (/-dependence in Sf" 1 , consider its 
derivative with respect to q: 

dSr(A d ) , ( k°, \ . . 

_^LJ = 4( J V-l)(2.r) tr ( M ^j. (2.23) 

The integral is most easily done by integrating first over the spatial momenta: 
+oc r rl 3 k ( k° \ + °° 

(V " L) T „£. I wf lwf+p) = ^ 

This sum, while formally divergent, is interpreted using zeta function regularization 



15 1 . The zeta function ((p, q) is defined as 



Hence 



Using 



{n + q) p 



tr ( 77l3 =* ) = -(Ytr L) ttT 3 (C(-2, q) - C(-2, 1 - q)) . (2.26) 



Xk%) 2 + k 2 



12 dq 

integration of ( |2.23| ) gives 

S! u (A d ) = V tr ^— (N-l) j dz q 2 (l - q) 2 . (2.28) 

From the nature of the sum in ( 2.24 ), Sf u is periodic in q, and is invariant under 
shifts of q — > q + I for any integer I. Thus in ( p. 28 ) q is defined modulo one. Also, a 
g-independent constant in (|2.28|) was dropped; this constant is just the free energy 
of an ideal gas of iV 2 — 1 gluons at a temperature T. 

As promised, the classical degeneracy in q is lifted by quantum effects: the minima 
of the theory are now at integral values of q = j, where trQ = exp(2irij / N) . 
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In (|2.28| ) the length in the z direction, L, is replaced by the integral over z. Of 
course for a constant field this substitution doesn't matter, but consider the interface 
problem. Introducing the dimensionless coordinate z', 




z' = J-gTz, (2.29) 
the sum of the classical and quantum actions becomes 



2 



S d {A d ) + ST{A d ) = V tr l) - [dz> [(£) +q 2 {l-qf\ . (2.30) 



f%N g J \\dz 

We can view minimization of this effective action as a problem in mechanics, with 
the coordinate z' as the "time". The classical action contributes the kinetic energy, 
(dq/dz') 2 , while the quantum action produces a standard double well potential, q 2 (l — 
q) 2 . For any solution to the equations of motion the energy, £ = (dq/dz') 2 — q 2 (l — q) 2 
is conserved, d£/dz' = 0. With the Z(N) interface we want a solution which obeys 



the boundary conditions q(0) = 0, q(L') = 1, as V = (^jN/3)gTL — > oo |5]||. By 
the boundary conditions the instanton has zero energy, £ = 0. Consequently, for the 
Z(N) instanton 




2 f 1 dq q(l-q) = \ . (2.31) 
Jo 3 



If the total action of the interface is the transverse volume Vj r times the interface 
tension a, then at leading order || 

a = T - . (2.32) 

3V3N g 

Two assumptions must be justified. The first is why Sf u can be computed for a 
constant field q, and then applied to the Z(N) instanton, where q is clearly a function 
of z. The reason can be seen from the definition of z' in ( |2.29| ). As there is no length 
scale in the rescaled action, in terms of z' the instanton's width is of order one. For 
the original coordinate, z, this implies that the Z(N) instanton is "fat", with a width 
of order 1 / (gT) . In weak coupling this is much larger than the natural length scale in 
a gas of massless, nearly ideal gluons, which is 1/T. Thus while the Z(N) instanton 
field is large in mag nitude, A d ~ T/g, it varies slowly in space, and at leading order 
this variation can be neglected. Corrections to this approximation do enter beyond 
leading order. 
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The second assumption is whether in the space of SU(N) gauge fields the path 
chosen is of minimal action. The Z(N) instanton interpolates between trQ — 1 and 
tr Q = exp(2iri / N) . By a global gauge rotation Q can be chosen as a diagonal matrix, 
involving the (N — 1) diagonal generators of SU(N) , the Vs. The quantum action 
for a general (constant) field can be computed directly [§]. We chose the simplest 
path possible — straight along the direction — but it is not obvious that other 
paths, which wander off into the direction of the other (N — 2) ij's, might not have 
lower action. (Our path is at least a local minimum.) For SU(2) there is only one 
path possible. In appendix A we show that for N = 3 and N = oo, the path along 
tjv is minimal. On this basis we conjecture that this remains true for all 4 < N < oo. 



III. General analysis of the Z(N) interface 



The partition function of an SU(N) gauge theory is 

[dA„{x)\ [d V (x)\ e ^'(A)- SS f { A, v) ? (3 y 



where S cl and S 9 ^ are the classical and gauge fixing actions of (|2.14j ) and (|2.16|) . The 
gauge field = + A q ^\ for the time being the choice of A^ is left open. The 
coordinate of four dimensional spacetime is x = (r, Xt, z) 

In order to reduce this four dimensional theory to an effective theory in one di- 
mension we introduce the field 

q(z) = -L J'dr J d 2 x t g <(*) . (3.2) 

This choice is obviously motivated by the definition of Aq\ unlike the scalar field q of 
( p.2[ ), q is a matrix valued field in the adjoint representation of SU(N) . The only 
subtlety in the introduction of q(z) is our insistence on defining it not just as the 
average of Aq over time, but over the spatially transverse directions as well. The 
reason for this will become apparent shortly. 

The field q(z) is introduced into the functional integral by a delta function con- 
straint: 

[dA,{x)\ [d V (x)} [d\(z)}[dq(z)} e -s^)-^(^)-s-(A,q,A) ; (3 3) 



S con (X,q,A) = J dz 2i tr (\(z) (q(z) - ^- dr Jd 



^A^x) 



(3.4) 
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The constraint field X(z) is introduced to enforce the definition of q(z) in ( |3.2| ), and 
so is also an adjoint matrix. With the overall factor of i in S con , the original contour 
of integration for the constraint field A (z) is along the real axis. 

As an aside, note that q(z) does not transform in a simple fashion under gauge 
transformations. For an infinitesimal gauge transformation ou, where — ► d^uj — 

q(z) -> i- j'dr J cPx t g (-ig)[A (x),u(x)} , (3-5) 

assuming that a; (a?) = u}(r,x t ,z) is periodic in r. Nevertheless, no further terms 
besides those in (|3.3|) are required in the measure of the functional integral: gauge 
fixing for takes care of that. 

The effective action S e ^(q), is defined as the integral over all fields, excepting 
qO): 

Z = J [dq(z)\ e - 5e//(q) . (3.6) 

How the effective action is computed in practice depends upon the problem at hand. 
But by introducing q(z) as an extra field into the functional integral, clearly there is 
no confusion possible about double counting degrees of freedom. 

The calculation of S e ^(q) is straightforward at one loop order J7|, and so our 
discussion is brief. The field q(z) is held fixed, while the gauge field and the constraint 
field A are seperated into classical plus quantum terms: A^ = A C J + A q ^ and A = 
A d + \ qu . Expanding the action of ( |3~3|) to quadratic order, 

S C \A) + S 9f (A, v ) + 5 con (A,q,A) ^ S + S 1 + S 2 . (3.7) 

The leading term is the sum of the classical action plus the constraint: 

S = S d {A d )+S con {X d ,q,A d ) . (3.8) 



The linear terms determine the equations of motion: 

S, = S con (X qu , q, A d ) + 2iJ d 4 xtr(\ d 5 v o + iD d G d u ) A q ^ . (3.9) 

The constraint term, S con (X qu , q, A d ), is linear in X qu , and so determines A cl (x). The 
obvious choice is to take A cl (x) to be a function only of z, 

?ttT 

A d {z) = — q(z) . (3.10) 
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The constraint term in Sq then vanishes, S con (X cl , q, A cl ) = 0. 

The second term in Si modifies the equations of motion for the gauge field: the 
constraint term acts as a source for the gauge field, proportional to X cl . As A cl is 
determined, let 

X d (z) = -iDfG%{z) , (3.11) 

from which S\ = 0. (As A cl (x) is assumed to depend only upon z, the other com- 
ponents of the equations of motion for the gauge field are automatically satisfied.) 
At the stationary point X cl is purely imaginary, while the fluctuations X qu remain 
real. This shift in the contour for the constraint field A — by an imaginary amount, 
keeping it parallel to the real axis — is standard. 

The quadratic terms in the action are 

S 2 = - [ d'x{2t) fr(A*»(z) Af{x)) 



2uTV f 



tr 



+ Jd'xJ d'y tr (A*tfjT x A? + V (-Dl) v ) . (3.12) 

The first term in S 2 > involving X qu , is special to the constraint action. Since A 9 " only 
enters into S 2 linearly, it can be integrated out. This introduces a constraint for the 
integration over Aq U (x): 

[ dr fd 2 x t A^{r,x t ,z) = 0. (3.13) 

Vtr JO J 

This constraint is completely innocuous. Remember that for the two spatial directions 
x t , each length L t is taken to infinity. Integration is over all modes of X QU (z) which 
obey ( |3.13| ), but the only modes which don't are those constant both in r and x% - 



in momentum space, modes with k° = k t = 0. If the length in the x t directions are 
infinite, the corresponding momenta kt take on all continuous values, and those with 
k t = have zero measure, and can be ignored. 

Being able to drop the constraint is an important point. Suppose the constraint 
field is defined not as an integral over r and xt, but just as an integral over r: then 
q is a field in three, instead of one, dimension. Going through the same proceedure 
as above, at quadratic order integration is over all fields constant in time. But these 
modes have nonzero measure in the functional integral. At finite temperature, the 
momentum A; = 2irnT for integral n, and the constant modes, with n = 0, are of 
countable extent. 
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For situations in which the constraint doesn't matter — that is, where the effective 
fields are of zero measure in the space of the original fields — the constraint methods 
give the same result as for the usual effective potential J7|. The constraint field X d 
plays the role of the external source, while the exchange of A cl for q mimics precisely 
the process of Legendre transformation. The effective potential is the energy of the 
vacuum in the presence of the external source, and so is properly minimized. 

The meaning of the effective action when the effective fields are of nonzero measure 
in the space of the original fields is unclear. We dwell on this point because of such 
an analysis by Oleszczuk and Polonyi [It|, who introduce an effective field in three 



dimensions, as the integral of Aq with respect to time. Integrating out modes with 
A; 7^ 0, they find a potential different from that of ( |2.28|) . The potential in (|2.28 



is a 



constant times q 2 (l — q) 2 = q 2 — 2q 3 + q 4 . The potential for the three dimensional field 
of JTOl is just q 2 + q . the term — 2q 3 is missing, as that arises from the k° = mode of 



the integral. More generally, with an effective three dimensional theory, the manifest 
Z(N) symmetry of the original theory is broken by the separation into modes with 
zero and nonzero k°. The loss of the Z(N) symmetry seems a grievous price to pay. 

Returning to ( p,12|) , the remaining terms are standard in a background field ex- 



pansion. In background field gauge, Q2.16|) , the inverse gauge field propagator is 

yd I 

fi — V 1 "J Y~" liVI J 1 



Ajr 1 = -Dl <$"" + DiD% - -D d D d + ig[Gl 



= ~D d ff" + f 1 - - J DfD? + 2ig[G%,} . (3.14) 
After integrating out X qu , A qu , and r), 

S eff (ci) = S c \A d ) + Sr(A d ) , (3.15) 
where A cl is related by q(z) by ( |3.10| ), and 

ST(A*) = \tr log (AjT 1 ) - tr log (-D 2 d ) . (3.16) 

For the Z(N) interface let 

q(z) = q(z) t N , (3.17) 

and then repeat the analysis of sec. II. Taking the field of the Z(N) instanton as 
slowly varying, to leading order in g 2 the effective action for q(z), S e ^(q), is given 
by (pi. 
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At next to leading order corrections arise from two sources. Viewing the action of 
( p.30| ) as a type of quantum mechanics, these terms can be understood as corrections 
to the potential and kinetic terms. As the potential was first generated by the free 
energy at one loop order, so corrections to this potential are produced by the free 
energy at two loop order. These two loop effects are g 2 times those at one loop 
order. Like the calculation at leading order in sec. II, for the two loop potential the 
background field can be taken as constant. Secondly, at next to leading order it is 
necessary to account for the spatial variation of the Z(N) instanton. For this only 
the free energy at one loop order is required, expanding to leading order in (dq/dz) 2 ; 
thus these terms correct the kinetic term in the effective action. With our definition 
of q, the classical action is proportional to 1/g 2 times (dq/dz) 2 , (|2.15|) ; the terms from 
the free energy are of order one times (dq/dz) 2 , and so are smaller by g 2 . 

Ultimately, corrections are small because the Z(N) instanton is fat: the ratio of its 
size to the thermal wavelength is of order 1/g. Hence an expansion in the derivatives 
of the instanton field is automatically an expansion in g 2 . This is what makes the 
problem tractable. 

Both of these effects are due to the effects of fluctuations in four dimensions as 
they generate the effective, one dimensional action S e ^(q). The functional integral 
over q in (|3]6|), however, is still treated classically. When do fluctuations in q(z) 
enter? 

Remember that the quantity of physical interest || is the partition function, Z, 
for a system with the appropriate boundary conditions to enforce a domain wall in 
the spacetime tube: 

Z = ce"* . (3.18) 

While we have concentrated on the interface tension, a, of course the prefactor "c" 
is also of significance. 

Integration over fluctuations in four dimensions generate the effective theory in one 
dimension, S e ^ , and so determine a. Fluctuations in S e ^ , though, do not contribute 
to a, only to the prefactor. This is simply because S e ^ itself is proportional to V tr , 
and so the integral over the effective, one dimensional fields cannot generate a constant 
times Vtr in the exponent, but merely powers of Vt r in the prefactor ||. The prefactor 
for a Z(N) domain wall is given by the integral over S e ^ at one loop order; however, 
it is necessary to compute S e ^ for a general path in group space, (|a. 1| ) , instead of 
the "classical" path of (|3.17|) . This we defer. 
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Consequently, we confess that the machinery of the effective action S e ^(q) devel- 
oped in this section is not essential for what follows. We discussed it at such length 
in order to ensure that there are no problems of principle, and to emphasize the 
generality of the method. 

IV. Z(N) interface at next to leading order: Feynman gauge 

In this section we compute the leading corrections to the interface tension in 
background field Feynman gauge, £ = 1; in the next section, for arbitrary £. While 
the methods are the same for all £, technically the calculations are simpler in Feynman 
gauge, and so in this section we discuss our methods in some detail. In sec. V the 
results are merely summarized, in order to emphasize the physical interpretation of 
the gauge dependence which arises in the effective action. 

As discussed following (|3.17|) , at next to leading order there are two pieces needed 
for the interface tension. The first is the effective potential in a constant background 
Aq field to two loop order. This was calculated for 577(3) by Belyeav and Eletsky ]ll| 
and by Enqvist and Kajantie ||12|| . We have independently computed the potential 
for the field of ( |2.2|) at arbitrary N, but given previous calculations, are content here 
to just establish the iV-dependence at two loop order. 

At one loop order the N dependence of Sl u (A cl ) is obvious. The background field 
enters only through adjoint covariant derivatives; from ( |2.11| ), only has nontriv- 
ial commutators with the generators t% y Thus the only g-dependence is from the 
free energy of the 2(iV — 1) fields for these ladder operators, and at one loop order 
Sl u (A d ) ~ (N-l), (gggp . 

The diagrams which enter at two loop order involve either two three-gluon vertices 
or one four-gluon vertex. Both types of diagrams involve a product of structure 
constants. The only diagrams that depend nontrivially upon the background field 
are those in which two lines are along the ladder operators t% r Denoting the q- 
dependent terms in the free energy at two loop order as Sf 1 (A cl ), after writing each 
structure constant as a trace 

N-i 2 2 

sr{A d ) ~ EE Hm,i,%])) = (N-i) E^^^-D'^-dO) • 

j=l a a 

(4.1) 

The sum is over all generators t a with nonzero trace. The last expression follows by 
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as j — (N — 1). The commutator for t^ ( - Ar _ 1 - ) is: 



V N,(N-1)' °N,(N-1)1 ~ 2 



(4.2) 



noting that from the form of t jy in (|2.3| ), each value of "j" contributes equally to 

1]). Thus the complete sum over "j" is (N — 1) times that for any single term, such 

± 

N,(N-l) 

( l ) 

... -1 

v ... l) 

All elements not indicated vanish. From (|4.2j ) the only terms which contribute to 
( |4. 1[ ) are if t a is one of two diagonal generators, or t^-i, ( |2.3j ) and (|2.5|). It is 
then easy to show that the sum over "a" in ( |4.1| ) is proportional to N, so that in all 
St u (A d ) ~ N(N-1). 

Knowing the iV-dependence, the general result can be read off from that for N = 3 
|I~2|| . For the constant field of ( |2.2| ), the sum of the the free energies at one loop 



order, Si U (A cl ) fl2.28Q , and at two loop order in Feynman gauge, S^ =l (A cl ), is 

Sl u (A cl ) + 5 2 ^ =1 (A d ) 

= V tr ^(N-l) J dz [q\l- q f + (0)(3? 2 (l-# " 2g(l-g))) . 

(4.3) 

The second piece required for the interface tension arises from the free energy at 
one loop order for a background field which varies in z. This corrects the kinetic 
term in the effective action, as a function of the background field, q, times [dq/ dz) 2 : 
once one factor of (dq/dz) 2 is extracted, the remaining factors of q can be taken as 
constant. At one loop order the quantum action Sf u (A cl ) is given by ( |3.16|) , with the 
inverse gluon propagator of Q3.14j ), setting £ = 1 in Feynman gauge. To calculate 
this some, although not all, of the tricks of sec. II can be used. For example, from 
( P-17| )-( p.l9| ), for a constant field the one loop quantum action is independent of £. 
This is no longer true for a spatially varying field. 

The ladder basis of sec. II can be used to simplify the color algebra. Defining 
Gq z = —G+Q = dq/dz, else zero, and D\ = {D^) 2 + d 2 , with Dq as in ( |2.20|) , for 
arbitrary fields q(z) the one loop quantum action reduces to 

Slt=M d ) = (N-l) {trlog(-D 2 + 5, v + ^TiG+ v ) - 2 tr log (- D 2 + )) . (4.4) 

One kinetic term arises by expanding to quadratic order in 

S^M") - - 16^ (N - 1) (g 2 ,r (po^) • (4.5) 
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with k° + as in ( ggg ). 

Having extracted this, the remaining kinetic term arises from the expansion of 

Slt i= M d ) = 2(N-1) (l - trM-^). (4.6) 

The factor of 1 — e/2 appears because the theory is regularized in 4 — e dimensions, 
and in covariant gauges the number of gluons equals the dimensionality. 

Calculating the momentum dependence for such a one loop action is a standard 
problem: see, for example, the treatment of Iliopoulos, Itzykson, and Martin ||17|| . The 
computations of [fTTH , however, are in coordinate space, which for most problems is 
rather awkward. Instead, it is much simpler to perform the calculations in momentum 
space. We have checked that for the problem at hand, as well as for the scalar example 
treated in |T7|], the results agree. 

To work in momentum space, let q — > q + 5q in ( |4.6| ), and expand to quadratic 
order in 5q. The idea is to isolate the momentum dependence in the fluctuation, 5q: 

trlog(-Dl) ~ 8tt 2 T 2 Q tr (-^ (8q)^j + ^(^F^^F^)) ■ (4-7) 

The momentum dependence only arises through the second term on the right hand 
side, since the first term is a type of tadpole, independent of the momentum flowing 
through 5q. The field q varies only in z, so its momentum is purely spatial. If (0, p) 
is the external momenta, then, and (k°, k) the loop momenta, ( |4.7|) becomes 

8tt 2 T 2 5q(p)5q(-p) tr ( S£ . (4.8) 

\((k° + ) 2 + k 2 )((k + ) 2 + (p-k) 2 )) 

In this form it is trivial to expand to order p 2 . Trading L (5q (p 2 ) 5q) for / dz (dq/dz) 2 , 
the second kinetic term is 



dz) V 2J \((k° + ) 2 + k 2 f 3-e ((k%) 2 + k 2 y) ' 

(4.9) 

At nonzero temperature dimensional continuation is carried out by changing the 
number of spatial dimensions to 3 — e, which produces the factor of 4/ (3 — e) above. 

The integrals required are given in appendix B, (0), (|7J), and (|E?D . Without 
worrying about their detailed form, one feature is evident. Each integral is loga- 
rithmically divergent in four dimensions, so in 4 — e dimensions, there are poles in 
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1/e. These are the standard terms which produce the renormalization of the coupling 
constant at one loop order. Thus it is instructive to combine the results for (|4.5|) and 
( |4.9|) with the classical action of (|2.15 ) to find 

S d {A d ) + SH^ x {A d ) + Slu=M d ) 

4tc 2 T 2 . . r ( dq\ 2 ( 11 g 2 N ( , , , , 1 \\ 
= V tr -——(N-i) dz [-?-) l + _|_ 



g 2 (T)N v W ^cfe/ V 3 16tt 2 V 11// ' 

(4.10) 

ip{q) = d(logT(q))/dq is the digamma function. The prefactor includes the running 
coupling constant at a temperature T, g 2 (T), which is related to the bare coupling 
constant g 2 as 

with // the renormalization mass scale. The relationship between the bare and renor- 
malized coupling constants in ( |4.11| ) is arbitrary up to a constant; our choice is similar 
but not identical to the modified minimal subtraction scheme, and is convenient at 
nonzero temperature. At high temperature, Q4.ll ) exhibits the standard logarithmic 
fall off of the running coupling constant, g 2 (T), with the coefficient of lliV/3 appro- 
priate for the /3-function of an SU(N) gauge theory at one loop order. Notice that 
the g-dependence in (|4.10|) , through the digamma functions of q and 1 — q, enters 
with precisely the same coefficient as for the /3-function, lliV/3. 

The effective action which governs the Z(N) instanton at next to leading order 
is the sum of (|4.3|) and (|4.11|) . The action is determined by the properties of the 



solution at leading order. As discussed following ( |2.30| ), the Z(N) instanton has zero 
energy, and so 

2 

(V>(g) + il>{±-q)) 




V(l _ (j]) .995018... ; (4.12) 

the value of the integral was determined by numerical integration. Hence at next to 
leading order, the interface tension a is 

ot = « N „-J* 2 £ (l - (15.2785...)^) , (4.13) 



3^3N g(T) \ v •" 16vr 2 
which is the result quoted in ( |1 . 1|) . We have taken the liberty of writing the corrections 
as proportional not just to the bare coupling constant, g 2 , but to the running coupling 
constant, g 2 (T). 
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V. Z(N) interface at next to leading order: general gauges 



The calculation of the interface tension in an arbitrary background field gauge is 
similar to that for Feynman gauge. Nevertheless, the calculation illuminates some 
features which are missed by working at fixed £. 

Including the terms at both one and two loop order, the free energy in the constant 
background field of (|2.2|) is 

Sr(A d ) + S%(A d ) = Vtr ^^(N-l) J dz (g 2 (l-g) 2 

+ (0) ((7 - 40 g 2 (l - qf - (3 - £) q(l - qj)\ . (5.1) 



The two loop potential for general £ was computed first by Enqvist and Kajantie W2 
we agree with their result when £ = 1, but not for £ 7^ 1. 

The potential changes in a rather dramatic fashion in going from one to two loop 
order. At one loop order the potential is just a standard double well, q 2 (l — q) 2 . At 
two loop order the potential becomes ^-dependent. This includes a correction to the 
coefficient of the double well potential, as well as a new term, proportional to g(l — q). 
This new term is peculiar, for it controls the behavior of the potential for small q. If 
£ < 3, the stable minima are not at q = and q = 1, but at go ~ (3 — £)g 2 and 1 — go- 
On the other hand, if £ > 3, the stable minima remain g = and g = 1. 

Such a nonzero value of the stable minima would have profound consequences for 
a gauge theory at high temperature |9|-|T^|. For N > 3, the trace of the Wilson line in 
the fundamental representation is a complex number. Under charge conjugation (C) 
or time reversal (T) transformations, the trace of the Wilson line goes into its complex 
conjugate (up to global Z(N) transformations). Thus if the stable vacuum indeed 
has g 7^ (modulo 1), then the vacuum spontaneously breaks C and T symmetries, 
conserving CT. While conceivable, it is unexpected to find C symmetry breaking 
arising spontaneously in a pure gauge theory. Of course a physical phenomenon 
cannot depend upon the choice of the gauge fixing parameter, while go changes with 

£• 

Leaving these questions aside for the moment, in a general background gauge the 
kinetic terms in the effective action are, to one loop order, 

s c \a c1 ) + sr a ,M d ) + s lU Ad ) 
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4vr 2 T 2 . . r , /dg\ 2 / 11 o 2 iV / . . , . 7-6£\\ 

(5.2) 

The running coupling constant, g 2 {T), remains as in ( f4.11|) . 

In ( |5.2|) the only £ dependence is an overall constant, proportional to 7 — 6£, and 
is independent of the background field q. As discussed following ( |b.7| ) in appendix B, 
for each of the integrals which contribute to the kinetic term, the coefficient of the 
pole in 1/e is always the same as for the digamma functions of q and 1 — q, which is 
how the g-dependence arises. As is customary in background field calculations [[L6|, 
the poles in 1/e generate the /3-function at one loop order, and is independent of £. 
Thus if the coefficient of the digamma functions is the same as for 1/e, at one loop 
order it also must be independent of £. The only remaining ^-dependence possible is 
as a constant, which does appear. 

The effective action for arbitrary £ is the sum of ( |5.1|) and (|5.2j ) . It is easy to show 
that while each term depends individually upon £, for any solution with zero energy, 
£ = 0, the ^-dependence cancels in the sum. As the Z(N) instanton has zero energy, 
the value of the interface tension for £ ^ 1 is equal to that for £ = 1, (|4.13|) . 



The cancellation of ^-dependence is a necessary check on the consistency of our 
method, but by itself is rather unsatisfactory. To understand this better we follow 
Belyaev [|14J]. In SU{2) Belyaev showed that the apparent ^-dependence in the po- 
tential for a constant Aq field can be understood as a renormalization of the Wilson 
line. We now generalize his results to arbitrary SU(N) , and show that they explain 
the ^-dependence of both the potential and kinetic terms in the effective action. 

The point is that while the vacuum expectation value of the trace of the Wilson 
line is a gauge invariant quantity, the fields which we have been using to parametrize 
the Wilson line — A$, and so q — are not. For instance, from (|3.5|) , our effective field 
transforms in a nonlocal manner under infinitesimal gauge transformations. At tree 
level this doesn't matter, but it does at one loop order and beyond, as the Wilson 
line undergoes both infinite and finite renormalizations. To compute these, define 

A (x) = A$ + Al u (x), (5.3) 



with Aq related to q as in (2.2). In general the Wilson line is a function of the 



spatial position, and so we consider the trace of the Wilson line, averaged over space. 



20 



Expanding in powers of A qu , 

1 IT I T ^ = tr + fi i + ^ + • • • • (5.4) 

The first term on the right hand side, tr Q(A cl ), is the value in the classical background 
field, (|2.4D . The term linear in A qu , can be written as 

fil = 19 I T l dT j^^(Af{T,x t ,z)n{A d )) . (5.5) 

Due to the constraint imposed upon the quantum fluctuations in ( |3.13| ), this term 
vanishes. 

The term quadratic in quantum fluctuations is nontrivial. Consider first its value 
in zero background field, A d = 0: 

g\N* - 1) r d*k 1 

n2 = — 2 — n (27F2F- (5 - 6) 

This is a standard renormalization of the Wilson line, a type of "wave function" 
renormalization, proportional to (3, the length in euclidean time. In dimensional 
regularization this vanishes identically. 

However these infinite terms are regularized, they are independent of the back- 
ground field. In addition, there are finite terms which depend upon A cl . Using the 
commutation relations of Q2.11D , for an arbitrary constant u y" 

e ytN tN tj e~ ytN = e ±ycN t%j ■ (5-7) 

Using this relation, and remembering the path ordering required for the Wilson line, 
the terms dependent on the background field are 

2 N-l „a fT i 

Q 2 = - 9 - £ I drj o dr> (AZ + (t) Aiy(r') e~ 2 ^-^ tr (t% Aj Q(A d 

+ KA T ) A oA T ') e2mq{T ~ T ' ] tr (Vh Q ( Ad ))) ■ ( 5 - 8 ) 

Aq" ± is the component of Aq U in the direction of t%j. The r integrals are evaluated 
using the background field propagator, as in ( |2.18| ). After doing the color trace, in 
momentum space 

Q 2 = ig2 ^ N ~ l) e^' N (l - e—) tr (±. A 00 (k° + ,k)) . (5.9) 
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A o(/c°, k) is the usual covariant gauge propagator, 

Aoo(*°, *) = - (1 - J^f^y • (5.10) 

with the only dependence on the background field in (|5.9|) through the shifted mo- 
mentum k° + , fl2~2Tl) . 

How then to interpret these corrections at one loop order to the classical value 
of the Wilson line? By the constraint imposed upon A qu , the linear term in (5.5) 
vanishes, and so they cannot absorbed in A qu . We then introduce "renormalized" 
fields for Aq and q as 

A r en = \ T q ren t N , q ren = q + 6q, (5.11) 
gcN 

and require that the spatial average of the vacuum expectation value of the Wilson 
line be given by A ren : 

I IF I T ( tr n ^ = tr n ( Aren ^ ■ ( 5 - 12 ) 

Expanding the right hand side to linear order in Sq, 5q is proportional to ^2- The 
integrals required for ( |5.9|) are (|b.9|) and (|b,10| ) of appendix B, and give 

Q 2 N ( 1\ 

q ren = q + (3 - (q ~ =: J ■ (5.13) 



167T 2 

This relation is valid for < q < 1. For SU(2), it agrees with the result of Belyaev 



| 14| . Note that the renormalization from q to q ren is entirely a matter of a finite shift. 



The gauge dependent actions, (|5.1| ) and (|5.2f) , can be trivially rewritten in terms 
of the renormalized fields. The potential term becomes 

47r 2 T 4 r ( o 2 N\ 

Sr(A-) + ST(A-) = V tr (N-l) J dz (l - 5 f^J q 2 (l - qf , 

(5.14) 

while for the kinetic terms, 

S c \A ren ) + SZ(A ren ) + S^(A ren ) 

(5.15) 
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With the effective action of ( |5.14| ) and ( |5.15| ), the corrections to the interface are of 



course unchanged, equal to (|4.11| ). 

Once written in terms of q ren , all of the gauge dependence found previously in 
the kinetic and potential terms cancels. Further, the effect of two loop terms in the 
renormalized potential is just to change the coefficient of the one loop terms: the 
new terms found previously at two loop order, proportional to q(l — q), have all been 
absorbed by q ren . 

Hence the apparent instability of the perturbative vacuum at two loop order 
merely results from a classical parametrization of the renormalized Wilson line. After 
correcting for loop effects, the stable vacuum is the trivial one (plus Z(N) transforms 
thereof) and is C symmetric. 

It seems unlikely that the cancellations found at two loop order are mere coinci- 
dence. We conclude with a conjecture: that the stable vacuum of hot gauge theories - 
both with and without fermions — is symmetric under charge conjugation to arbitrary 
loop order. 

The research of A.G. and R.D.P. was supported in part by the U.S. Department 
of Energy under contract DE-AC02-76-CH0016. 

Appendix A: Proof of minimal action for N = 3, oo 

In this appendix we prove that for N = 3 and oo, the path chosen for the Z(N) 
instanton is of minimal action. 

By a global gauge rotation a constant background field A d can be chosen to be a 
diagonal matrix. Thus the most general constant field for the A cl of (|3.10|) is 



The ti are the N — 1 diagonal generators of SU(N) , as in (|2.3j ) and (|2.5[) . For clarity, 
in ( |a. 1| ) the field q is relabeled as qw- 

Promoting each qi to be a function of z, for this ansatz the classical action becomes 
a sum over N — 2 independent kinetic terms 



N 



q = £ Qi U ■ 



(a.l) 



i=2 




(a.2) 
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For the Z(N) interface the boundary conditions required are 

g»(0) = gi (L) = , i = 2...(iV-l); q N (0) = , q N (L) = 1. (a.3) 

We wish to show that the path with q 2 (z) = . . .qN-i(z) = 0, and qn{z) = q(z) 
as before, is the path of minimal action. To demonstrate this we need the potential 
generated by fluctuations at one loop order. We do so in two special cases where the 
analysis is elementary. 

N=3: There are two independent fields, q2 and q 3 . The potential term is 

Sl u {A d ) = V tr ^- JdzV tot (q 2 ,q 3 ) , 

Vtot(q2,qs) = V(g 2 ) + V(q 2 /2 + q 3 ) + V(-q 2 /2 + q 3 ), (a.4) 

where 

V(q) = [q] 2 (l - [q]) 2 , [q] = \q\ modl . (a.5) 

Because the V(q) is a function of [q], the absolute value modulo one, it is not quite 
the simple polynomial form it first appears to be. (This restriction could be ignored 
before, since the Z(N) instanton only involves q : — > 1.) 

The effective action is the sum of (|a.3| ) and ( |a.4|) . It is not difficult to see why 
for the boundary conditions of ( |a.3|) , the path with q 2 {z) = is of minimal action. 
Considering the problem as classical mechanics in two dimensions, from energy con- 
servation the action for any solution to the equations of motion is proprotional to 

/ ds \IV to t{.q2,q?>) , (a.6) 

where ds is the arc length in the space of q 2 and q 3 . It can be shown that for any 
fixed value of q 3 , the potential is minimized for q 2 = 0: V to t(q2, ^3) > V tot (0, q 3 ). Given 
the boundary conditions, the path chosen is clearly of minimal length. As both the 
arc length and the potential are bounded by the path with q 2 = 0, by ( |a.6| ), so is the 
action. 

N = 00: The large iV limit is taken by holding g 2 = g 2 N fixed as iV — > 00; we work 
in weak coupling, for small g. Then the interface tension is of order N, a ~ N/g, 
while the Z(N) instanton remains fat, with a width of order l/(gT). 
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Besides q^, there are of order N q^s which can contribute. Assume first that there 
are a finite fraction of the q^s for which qi(z) 7^ 0. From ( |a.3|) , since each field has 
a kinetic term proportional to N, if order N fields contribute, the sum of the kinetic 
terms for all fields is of order N 2 . Similarly, the potential term is also of order iV 2 . 
Because both the kinetic and potential terms are positive definite, any solution has 
positive action, equal to a pure number times N 2 . This is N times the action of the 
path chosen and so is not minimal. 

Hence we can assume that at large N, only a finite number of the q^s contribute. 
For simplicity, assume that they are just qx and q^-i- The generators for t^, (|2.3|), 
and t JV— 1 5 (P-5|), simplify greatly at large N, reducing to just one diagonal element. 
(Essentially, the terms which enforce tracelessness can be ignored, as a correction in 
1/iV.) In this case it is easy to work out the full potential term, which is proportional 
to 

Vtot(qN-uqN) = (N-2)V(q N ) + (N-2)V(q N ^) + V(q N ^ - q N ) 

~ N(V(q N ) + V{q N -!)) . (a.7) 

The last term is the leading term at large N; like the kinetic energy, it is of order 
N. But notice that at large N, in the potential the coupling between q^ and qN-i 
has dropped out. Thus even though the potential for qN-i is nontrivial, because it is 
positive definite, the path of minimal action which satisfies gjv-i(O) = qN-i(L) = is 
just q N -i{z) = 0. 

This argument generalizes immediately to any finite number of gj's: the term of 
order iV in the potential is a sum over decoupled potentials, and the path of minimal 
action excites only one field, q^- 



Appendix B: Integrals in a constant A field 



In this appendix we catalog the integrals required for the computation of quantum 
actions in the constant background field of ( |2.2| ). 

The trace is defined as 

Remember that the system is of length L t in the x t directions, and L in the z direction. 
At a temperature T, (3 = 1/T, and the momentum k° = 2imT for integral n. The 
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volume of spacetime V = (3L 2 L. Dimensional continuation in 3 — e dimensions is 
used. 

As discussed following fl2.22|) , after using the ladder basis to reduce the color 
algebra, the propagators in the background field of ( |2.2| ) are given by their values 
in zero field, except for the replacement of k° — > k+ = 2nT(n + q). Following the 
example of (|2.22| )- (|2.28|) , the integrals are most easily computed by zeta function 
techniques First, the spatial integrals over d 3 ~ e k are done using the standard 

formulas of dimensional regularization. This leaves a sum over n, which is evaluated 
in terms of zeta functions. For instance, the simplest integral which arises in the 
expansion of the kinetic term at one loop order is 

1 V / e\ / , , „„x +2? 1 



^unrTPr - t^-( 1 + i) ^ ~ ,0 ^ T2 » „£rt^' (b - 2) 

Using the definition of the zeta function, (|2.25|) , and then expanding to order e, 

+oo 2 

E | T i 1+e = C(i + e,g) + C(i + e,i-g) (V>(<z) + , (b.3) 

with ip{q) = d(logT(q))/dq the digamma function. Hence in all 

1 V (2 



tr 



where 



Similarly, 



l = \+ V(V2) - log{*T 2 ) . (b.5) 



tr nuoVlw =7^(1- M«) + ^(1-?)) + 2) , (1..6) 
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These integrals are all those required to compute the one loop corrections to the 
kinetic term in see's IV and V. The factors of the renormalization group scale \i 
which enter into the running coupling constant, ( |4.11| ), arise by taking g 2 — > g 2 fi e in 
the classical action, and expanding in e. 

Notice that for the three integrals of (|b.4j) , ( |b.6j ), and (|b.7|) , the dependence on the 
background field q, through the digamma functions, always has the same coefficient 
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as the pole in 1/e. This is because in all of these integrals, q enters through the zeta 
functions £(1 + e, q) + C(l + e, 1 — q). Expanding about small e, ( |b.3|) , generates a 
common factor of 1/e and digamma functions in each case, so that up to their overall 
normalization, these integrals only differ by a constant. 

Other integrals required for the potential at two loop order, and the renormaliza- 
tion of the Wilson line at one loop order, include (|2.22|) , (|2.26| ), and the following. 
These are all finite integrals, so it is safe to set e = 0. 

tr MFTp = — (« 9 + e) • (b ' 8) 

1 VT ( 1\ 

tT fc° ((fc°)2 + fc2) = V ~ 2 ) ' (b ' 9) 

tr (MFW = " 17 I 9 " 2J ■ (b - 10) 

The last two integrals, ( |b.9|) and ( |b.lO| ), are singular for q — or 1, and as written 
are defined over < q < 1. Properly defined, they are discontinuous and vanish for 
q — or 1. 
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